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Abstract. We examine the gyration motion of a charged particle, viewed from a 
reference observer falling along the Z axis into a Schwarzschild black hole. It is assumed 
that the magnetic field is constant and uniform along the Z axis, and that the particle 
has a circular orbit in the X-Y plane far from the gravitational source. When the 
particle as well as the reference observer approaches the black hole, its orbit is disrupted 
by the tidal force. The final plunging velocity increases in the non-relativistic case, 
but decreases if the initial circular velocity exceeds a critical value of ~ 0.7c. This toy 
model suggests that disruption of a rapidly rotating star due to a velocity-dependent 
tidal force may be quite different from that of a non-relativistic star. The model also 
suggested that collapse of the orbit after the disruption is slow in general, so that the 
particle subsequently escapes outside the valid Fermi coordinates. 
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1. Introduction 

A star closely orbiting a black hole may be disrupted by a tidal force, most likely resulting 
in highly luminous flares. Observation of such flares can thus provide important clues 
to the physics of black holes (e.g, jlj). The high luminosity phase is of short duration, 
making detection rare, but recent X-ray observations 00] suggest a declining phase after 
the tidal disruption in a non-active galactic nucleus. A numerical simulation of tidal 
disruption and subsequent dynamical evolution has been carried out using smoothed 
particle hydrodynamics |I] , in which the typical velocity of the fluid elements was mildly 
relativistic, less than 0.1c. 

Recently, Chicone and Mashhoon j3] pointed out a very interesting phenomenon of 
tidal forces. When an ultra-relativistic velocity is considered, the tidal force exhibits an 
unexpected property contrary to non-relativistic dynamics. Suppose that an object of 
finite size falls towards a black hole along the Z axis. The object would be stretched 
along the Z axis and compressed along the X and Y axes under the Newtonian tidal 
force. The tidal force would therefore cause anisotropic acceleration and deceleration in 
the particle dynamics near a freely-falling observer. That is, from the viewpoint of the 
observer, out-going particles on the Z axis would accelerate, while out-going particles 
on the X or Y axes would decelerate. This is a non-relativistic expectation. Chicone 
and Mashhoon however showed that when the particle velocity exceeds a critical value 
cj\f2 ~ 0.7c, the particles would decelerate along the Z axis, but accelerate along the 
X or Y axes. Chicone, Mashhoon and Punsly [H] also discussed the case of a spinning 
particle orbiting a rotating black hole. In their one-dimensional model, the relativistic 
particles instantaneously escape from a valid region of the local observer frame. 

Here, we present another concrete model of the tidal effects on relativistic motion. 
We consider motion of a charged particle in a weak magnetic field, and Schwarzschild 
metric. Fermi normal coordinates are constructed around the reference observer falling 
into a black hole. The magnetic field is assumed to be uniform along the Z axis in this 
frame. Far from the gravitational source where the tidal force is negligible, the orbit 
is circular in the X-Y plane. As the reference frame approaches the black hole along 
the Z axis, the tidal force becomes important. But how does it disturb the circular 
motion? In particular, how does the relativistic rotating velocity affect the result? The 
spatial orbit of the particle is always bound, so that it is possible to pursue the long 
term evolution in the Fermi system. This toy problem may deepen our understanding 
of the tidal force. 

This paper is organized as follows. In section 2, the equation of motion of a charged 
particle in Fermi coordinates is discussed. In section 3, the tidal effect in the non- 
relativistic limit is discussed. In section 4, motion with relativistic velocity is considered, 
and numerical results are given. Section 5 is the discussion. 
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2. Basic equations 

We consider a reference observer whose world line is geodesic. Fermi normal coordinates, 
(cT, X k ) in a neighborhood of the observer can be constructed such that the world line of 
the observer is given by the observer's proper time T and X k = 0. The metric accurate 
to second order in X k takes the formj^j: 

<7oo = -1- F R Qi0j X i Xi + ■ ■ ■ , (1) 
g 0i = -~ F Ro jik XiX k + ---, (2) 

9ij — $ij ~ g F RikjiX k X l + • ■ ■ , (3) 

where F R a f3 1 s are components of the Riemann tensor along the observer. We 
consider the dynamics of a particle with mass m and charge q, labeled X h (T) or 
(X(T),Y(T), Z(T)) in the Fermi system. The equations of motion for the charged 
particle subject to the Lorentz force are jH] 

X* + ( F Ro l0j + 2 F R lkj0 X k + 2 F R okj0 X*X k 

+ F R ikjl X k Xi + 1 F Ro kjl XiX k Xi)xi = ^(A l - A°X*) , (4) 
where the overdots denote differentiation with respect to time cT, and 

r- = i - + 'R*,** + 1 ' w'*** + 5 W**, (5) 

A a = ^-F a(3 U . (6) 
m 

In eq.®, F al3 is the Faraday-Maxwell tensor and U a = (T, FX 1 ) is the four velocity 
of the charged particle. We assume for simplicity that the electric field is zero and that 
the magnetic field is constant along the Z direction in the reference observer system. 
For motion of a particle in a uniform magnetic field B, the electromagnetic acceleration 
(JHJ) is given by 

(A , A 1 , A 2 , A 3 ) = (0, ujTY, - uTX, 0) , (7) 

where uj = qB / (mc) is the cyclotron frequency. 

Using an explicit form of the Riemann tensor in Schwarzschild spacetime, eq.pj 
can be reduced to 

X + Kx(l- 2X 2 + -Y 2 - -Z 2 \ - -KX (hYY - 7ZZ) = +^;Y, (8) 



3 3 ) 3 V > cT 

Y + KY ( 1 — 2Y 2 + ^X 2 - ^Z 2 ^j - ^KY (5XX - 1ZZ) = ~^X, (9) 

3 3 



Z - 2KZ ( I - 2Z 2 + \x 2 + It 2 ) - -KZ ( XX - ) ')' ) = 0. ( I 0) 

where 



GM , . 

« = ^. (id 
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and r is the Schwarzschild radius, G the gravitational constant and M the mass. For 
a freely falling body with zero radial velocity at infinity, r can be expressed as a 
function of proper time T of the reference observer: r = (9GMT 2 /2) 1//3 , and therefore 
K = 2/(9c 2 T 2 ). We consider only the negative region of time T, in which T = — oo 
corresponds to r = +00, and T = corresponds to r = 0. However, when the spatial 
range is limited, time is also limited, e.g, when the spatial range is limited to r > r , 
time is limited to T < 



I2rl/{%GM). 

Equation (jlUj) shows that if the particle is initially moving on a Z = plane, (i.e. 
Z = and Z = 0), the particle always remains on the plane. From now on, we ignore 
motion in the Z direction. 



Non-relativistic limit 



In this section, we consider the non-relativistic limit of eqs.( 
spacetime. The relevant equations are 



and 



X + KX 
Y + KY -- 



+ -Y, 

c 

- u -x. 

c 



The analytic solutions of eqs. (|T2"j) and (fT3j) are expressed by Bessel J 
Ni/g functions of order 1/6: 



and assume flat 

(12) 
(13) 

and Neumann 



X 



Y 



2ixrj 



2nr] 



A cos 



V + X 



A sin 77 + x 



7T 

3 
3 



Ji/e(v) - 



- B sin yq + ip — 
B COS I 7] + ip — 



1/6 
7T 

3 

7T 

3 



,(14) 
,(15) 



where rj = lu\T\/2 and A,B,x,ip are integration constants. The meaning of these 
parameters becomes clear by taking the limit as infinity. For large |T|, we have 



X 
Y 



-A 
-A 



cos 



7T 



u\T\ + x + - ) - cosx 



sin uj\T\ 



X+ 3 



smx 



+ B 
+ B 



cos 



sin 



7T 



uj\T\ + ip + — ) + cos-^ 



7T 



uj\T\ + ip + — ) + sin?/' 



(16) 
(17) 



The orbit for large |T| represents a circle with radius (A 2 + B 2 — 2 AB cos(x — i ) )) 1 ^ 2 
and origin (Acosx + B cosip, A sin x + -Bsin^). Parameterization with A, B,x, *l> is 
awkward for large |T|, but useful for small |T| because two solutions, referred to as 
solution A (A 7^ 0, B = 0) and solution B (A = 0,B ^ 0), are remarkably different. 



Since Ji/ein) ~ rf^ and Ni/ e (n) 



rj 



-1/6 



for small rj, we have X, Y oc A x (77) 



2/3 



1,7 ocAx (?7)~ 1/3 in solution A, while X, F oc B x (r/) 1 / 3 , X,F oc 5 x (r^)" 2 / 3 
in solution B. Both solutions imply that the orbits converge at X = Y = 0, with 
the velocity diverging as |T| tends to zero. The velocity in solution B increases rapidly 
compared with that in solution A. General solutions are expressed by linear combination 
of the two solutions. 

In Figure 1, we demonstrate the trajectory of a charged particle with uo > 
experiencing a tidal force. The orbit at large \uT\ is circular with Larmor radius r c , 
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Figure 1. Orbits of charged particles subject to a tidal force are shown on the X-Y 
plane. Length is normalized with respect to the Larmor radius r c . The gyro-center 
at the initial gyro-motion is located at (r c /2,0) in the left panel and (3r c /2,0) in the 
right panel. 



while the final plunging orbit depends on the phase in the circular motion at a certain 
time. The figure shows seven orbits whose phase is given by 27m/7 (n — 0, • • • , 6) . The 
gyro-center of these orbits is fixed as (r c /2, 0) in the left panel, but is shifted to (3r c /2, 0) 
in the right panel. The tidal force causes the circular orbits to collapse, irrespective of 
the location of the gyro-center, and be compressed to the origin (0, 0) as T tends to 
zero. 




0.5 1 1.5 

-log (-C0T) 



0.5 1 1.5 

-log (-C0T) 



Figure 2. The amplification of velocity is shown as a function of time. The gyro- 
center of orbits is (r c /2,0) in the left panel, but (3r c /2,0) in the right panel. 

The velocity c(X 2 + Y 2 ) 1 / 2 is shown as a function of time in Figure 2. The value is 
normalized by the initial circular velocity, v = ur c . The velocity is constant during 
the circular motion, but significantly increases from a certain time T* in the final 
plunging phase. The transition time T* is ujT* —1, obtained by comparing the orders 
of magnitude of the Lorentz force qBujr c /c and the tidal force ~ GmM£/r 3 , where 
r(~ (GMT 2 ) 1 / 3 ) is distance from a gravitational source with mass M and the typical 
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distance £ from the reference observer is assumed to be of the same order as the Larmor 
radius, £ ~ r c . The tidal force dominates over the Lorentz force for small r, i.e, small 
|T|, and the condition of the tidal disruption can be approximated as \ujT\ < 1. 

After the transition time T*, the velocity increases. The curves in Figure 2 are 



of uoT — > 0. The non-relativistic treatment therefore becomes invalid, and the special 
relativistic effects should be taken in account in the final phase. 

4. Numerical results for relativistic case 

In this section, we consider the tidal effect on gyration of a relativistic charged particle. 
We neglect electro-magnetic radiation from the particle for simplicity. The tidal 
acceleration is described in Fermi normal coordinates constructed around a reference 
observer. The spatial range valid for the coordinate system is limited by radius 7£ -1 / 2 , 
where 1Z is a typical component of curvature tensor evaluated along the observer's world 
line. A simple criterion is 



Outside this region, the spatial metric in the X-Y plane changes sign. 

We assume that the center of the circular motion is located not far from the reference 
observer at the origin, i.e, X ~ r c , Y ~ r c . It is possible to realize \X\ ^> r c , \Y\ ^> r c , 
but the location of the reference observer is inadequate in such a case. In the non- 
relativistic approximation, circular velocity v in the absence of the tidal force satisfies 
the relation v = u>r c . Therefore, the location of the circular motion may be chosen to be 
close to the origin, as long as the velocity is small enough. Condition (j!8j) is therefore 
satisfied. Formally, the right hand side of eq. (jl8|) is infinity because the light velocity is 
assumed to be infinity in non-relativistic cases. 

In relativistic cases, the circular velocity satisfies v (1 — (v /c) 2 ) -1 / 2 = cur c . Therefore, 
X ~ r c , Y ~ r c , are not infinitesimal for v = 0(c) and a finite value of to. Far from the 
gravitational source, condition (fT%|) may be satisfied. As the particle falls, the curvature 
increases, and the valid spatial range accordingly decreases. As it is difficult to know 
a priori whether or not the particle always moves inside the valid spatial range, it is 
necessary to check condition (j!8|) in the evolution. To do this, we numerically integrate 
eqs.(jSJ) and Q with Z — 0, for relativistic orbiting particles. Orbits of relativistic 
particles are calculated from large |T|, provided condition (fTHJ) is satisfied. 

Numerical results are given for the motion with an initial circular velocity v = 0.1c. 
The trajectories are shown for the gyro-center (r c /2, 0) in the left panel and for (3r c /2, 0) 
in the right panel of Figure 3. Condition (fTHj) is no longer satisfied at the end of each 
curve, which corresponds to different time T. The final phase is therefore masked 
because the particle eventually escapes outside valid Fermi coordinates. This behavior 
seems to be general, because the particle motion is slow compared with the shrinkage 
of the valid spatial range, which is proportional to time |T|, as shown in condition (fTHJ) . 



only shown for values of uiT up to —10 3 , but actually diverge as 



|cuT| 2 / 3 in the limit 
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Figure 4. The same as in figure 2, but the initial circular velocity is 0.1c. 



The orbits are almost the same as those of the non-relativistic cases shown in Figure 1, 
the main difference being the termination points; this is due to the coordinate system. 
The relativistic effects on the orbits are not apparent. 

The locally measured value of velocity v at the particle position is given by 
v = c(— gijX % Xi / goo) 1 ! 2 , which can be reduced to c(X 2 + Y 2 ) 1 / 2 in flat spacetime. 
The time evolution of v normalized by the initial velocity v = O.fc is shown in Figure 
4. Like the non-relativistic cases, the velocity increases with ooT from the transition 
time ooT* — f. However, the final behavior, i.e, divergence of velocity seen in non- 
relativistic cases cannot be treated in this system, because the particle position is outside 
the valid coordinate system. The sudden drop in velocity in the right panel of Figure 
4 corresponds to the kink in the orbit shown in Figure 3. This is caused by the tidal 
acceleration almost reversing the direction of the particle's orbit. 

It is very interesting to consider the tidal effect on a relativistic particle with a much 
larger initial velocity. Is the particle accelerated up to the velocity of light? We show 
the results for initial circular velocity, v = 0.9c. The orbits are shown for gyro-center 
(r c /2, 0) in the left panel and for (3r c /2, 0) in the right panel of Figure 5. The spiral of 
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Figure 5. The same as in figure 1, but the initial circular velocity is 0.9c. 
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Figure 6. The same as in figure 2, but the initial circular velocity is 0.9c. 

the orbits is much tighter due to the larger circular velocity. In their final phase, the 
particles always escape from the valid coordinate system, as in Figure 3, preventing us 
from following the subsequent evolution. This termination may artificially cause these 
orbits to converge to (0.5r c , 0.2r c ) in the left panel and (1.5r c , 0.5r c ) in the right panel. 
Each position suggested from the figure is not a true convergent position because the 
end points do not correspond to the same coordinate time. Curves that stopped at early 
have larger (X 2 + Y 2 ) x l 2 due to condition (JTSJ). 

The time evolution of velocity v = c(— gjjX*X ? /g o) 1 ^ 2 is shown in Figure 6. It is 
normalized by the initial velocity v = 0.9c. Remarkably, the velocity never increases, 
and actually decreases. In this case, the tidal force causes deceleration, contrary to 
our intuition. The physical mechanism is not easily understood, but the mathematical 
origin is simple. The tidal acceleration terms coupled with velocity in eqs. (JHI) and © 
change sign for large velocities of X ~ Y ~ 0(1). Beyond a critical value, the tidal force 
becomes repulsive in the relativistic cases, whereas it is attractive in the non-relativistic 
cases. This property is analogous to that reported by Chicone and Mashhoon [2]. We 
have calculated evolution for various values of the initial circular velocity in order to 
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determine whether the tidal force causes acceleration or deceleration. It was found that 
the critical initial velocity is approximately ~ 0.7c, and that the velocity decreases above 
this value. 

5. Summary and discussion 

In this paper, we have examined the tidal effect on a charged particle orbiting in a 
uniform magnetic field. From this simple toy model, we found that the tidal force 
changes the property when velocity is relativistic, as shown by Chicone and Mashhoon 
jS]. In our model, the gyration velocity on the X-Y plane changes due to the tidal 
force coupled to the velocity, as the particle falls into a black hole along the Z axis. 
The velocity at the final plunging phase increases, when the initial velocity is less than 
a critical value. On the other hand, the velocity decreases, when it is larger than the 
critical value. The critical value was about 0.7c in our model, but depends on the 
process. 

As for speculative effects of the velocity-dependent tidal force in a more realistic 
astrophysical situation, suppose that a spinning neutron star falls into a black hole along 
the Z axis. The property of the tidal force depends on the rotating velocity v^, which 
corresponds to the orbiting velocity of the particle in our model. When exceeds a 
certain value, the tidal disruption may occur in a peculiar way, quite different from that 
in non-relativistic cases. The maximum rotation velocity of a neutron star is of the 
order ~ 0.5c at the equatorial plane (e.g, jH]). Therefore, it is very interesting to 
investigate the evolutionary behavior and the critical velocity in the hydrodynamical 
process with the velocity-dependent tidal force. 

As for the valid range of the coordinate system, the Fermi normal coordinates are set 
up around a geodesic reference observer, we used second-order expansion for the metric 
about the reference observer, and neglected higher order terms. The tidal force appears 
due to non-uniformities in the gravitational field, and is expressed by the deviation from 
the flat spacetime. The force becomes significant when the curvature tensors measured 
by the reference observer are large. At the same time, the higher order terms are no 
longer negligible. In other words, the valid range of the coordinate system shrinks as 
the gravitational source is approached. In order to study the final phase of the tidal 
disruption, it is important to include the higher order expansion terms. See e.g., [HUH] 
for the third order and fourth order expansion terms. 
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